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by 

Mary  Alice  Genau 


ABSTRACT:  Fourier  spectra  were  computed  on  the  IBM  7090  for  analog  tape 
recorded  pressure  pulses  of  underwater  explosions  fired  at  sea  in  February 
1965.  Depths  ranged  from  500  to  lU,000  feet;  charges  weighed  1  to  88 
pounds;  the  compositions  fired  were  T NT,  pentolite,  HBX-3,  and  Nitramex. 
Reduced  spectra  of  charges  weighing  up  to  57  pounds  agreed  with  previous 
results  fro®  1  and  10  pound  charges  at  the  same  depths.  Only  slight 
differences  due  to  composition  were  found. 
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SPECTRUM  AND  OCTAVE  BAND  ANALYSIS  OF 
PRESSURE  PULSES  FROM  DEEP  UNDERWATER  EXPLOSIONS 


1.  INTRODUCTION 

Weston  (reference  a}*  derived  an  analytical  approximation  of  a 
frequency  spectrum  analysis  which  is  applicable  to  the  pressure  pulse 
near  a  relatively  shallow  underwater  explosion.  In  order  to  determine 
the  effect  of  depth  upon  the  spectrum,  Christian  and  Blaik  (reference  b) 
performed  a  spectrum  analysis  in  19 6k  on  data  obtained  earlier  by  several 
different  investigators.  The  data  they  used,  which  was  all  recorded  near 
the  surface,  was  from  1-  to  10- lb  charges  fired  at  depths  between  7000 
feet  and  22,000  feet. 

To  extend  the  range  of  depths  to  shallower  conditions,  and  to 
include  larger  charges,  a  series  of  shots  was  fired  at  sea  in  1965.  The 
experimental  conditions  for  these  shots  are  given  in  Table  1;  for  compari¬ 
son,  the  conditions  for  the  earlier  experiments  are  also  shown.  This 
paper  reports  the  results  of  the  frequency  spectral  analyses  performed 
on  the  1965  shots. 


2.  EXPERIMENTAL  CONDITIONS 

The  experiments  were  carried  out  in  February  1965  about  200  miles 
east  of  Eleuthera,  in  water  about  18,000  feet  deep.  The  gage  for  recording 
the  pressure  pulses  was  suspended  from  the  USNS  GILLIS  at  a  depth  of  about 
200  feet.  The  charges  were  fired  at  depths  ranging  from  approximately 
500  to  lH,000  feet  directly  beneath  the  ship. 

One  Atlantic  Research  Corp.  Type  LC-3?  hydrophone  was  used  to  pick 
up  the  pressure  pulses.  The  same  gage  output  was  recorded  on  two  oscil¬ 
loscope  channels,  using  different  gains  and  sweep  speeds.  Almost  immediate 
observation  of  pressure-time  data  was  obtained  from  Polaroid  prints.  FM 
magnetic  tape  recordings  were  also  obtained  from  the  same  gage  at  60  inches 
per  second  on  the  FR-600  tape  recorder.  The  frequency  response  of  the 
EM  recorder  was  essentially  flat  from  0-20  kcs;  the  hydrophone,  however, 
gave  setie  low  frequency  distortion  (reference  c). 

Most  of  the  charges  fired  were  TNT  and  HBX-3;  1-lb,  8-11,  and  55-lb 
charges  were  used.  In  addition,  some  1-lb  pentolite  and  88-lb  Nitramex** 
charges  were  used.  Attempts  to  fire  1000-lb  TNT  charges  failed  (reference 
d).  All  charges  were  boostered  with  pentolite  which  in  turn  was  initiated 
with  hydrostatic  firing  devices  which  yere  set  for  nominal  depths  of  500, 
300,  1200,  3000,  4500,  7000,  10,000,  and  14,000  feet. 

*  References  are  listed  on  page  8. 

**  Manufactured  by  E.  I.  IXipont  de  Nemours  &  Co,  The  composition  is 
16*  TNT/ty  DNT/3C^  NaNO^/35^  NH4N03  plus  iron  and  phosphorus. 
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The  shot  data  are  given  in  Table  2. 


3-  DATA  ANALYSIS 

3«1  Method  of  Data  Reduction.  The  analog  tape  data  was  d"1  jiv-zed  by  the 
Mathematics  Department  of  the  David  Taylor  Model  Basin  on  their  Cotnpu*-  r 
Data  Format  Translator  (CDFT),  which  has  a  capability  of  '-'vopling  2000 
times  per  second.  Hie  analog  tapes  were  played  back  on  the  DTMB  Ampex 
FR-600  at  3-3 A  inches  per  second  and  sampled  and  digitized  1875  times 
per  second  on  the  CUT.  Since  the  data  was  recorded  at  60  incites  per 
second,  this  is  equivalent  to  sampling  at  33  microsecond  intervals  in 
real  time.  Hie  sampling  should  have  been  performed  at  25  microsecond 
intervals  in  order  to  observe  the  Jfyquist  criterion  of  sampling  20  kc  data. 
However,  sampling  at  the  next  lower  playback  speed  of  1-7/8  inches  per 
second  resulted  in  too  much  noise  and  so  was  not  done. 

A  fiducial  marker  one  millisecond  ahead  of  the  pressure  pulse  on 
the  analog  tape  initiated  digitizing  by  the  Compter  Data  Format  Translator. 
The  1-millisecond  baseline  was  digitized  so  that  the  average  value  of  the 
baseline  could  be  obtained.  Hv  sampling  rate  was  sufficient  to  average 
out  the  dominant  high  frequency  noise;  however,  the  baseline  was  too  short 
to  average  out  the  inherent  60  cycle  noise. 

Hie  number  of  times  a  pulse  was  sampled  ranged  from  about  500  samples 
for  the  deep  shots  to  5000  samples  for  the  shallow  shots.  The  digitized 
tapes  were  the  input  to  the  IHM  7090  computer  program,  NEWGRL,  described 
in  reference  (e). 

3*2  Types  of  Analysis  Performed.  Three  types  of  computation  were  carried 
out  on  the  IB4  7090  computer.  These  were: 

(1)  The  energy  spectral  density  E(f)*  of  the  positive  phase  of  the 
shock  wave  only  was  computed  as  described  in  reference  (e).  The  energy 
spectral  density  was  computed  in  increments  of  50  cps  from  50  cps  to 

l6  kes,  and  the  points  were  connected  by  straight  lines  by  the  CAL  COMP 
565  plotter . 

(2)  The  spectrum  was  similarly  computed  for  the  pressure  pulse 
through  several  pressure  oscillations  until  the  pressure  returned  to  the 
noi3e  level.  The  frequency  interval  at  which  the  spectrum  was  computed 
was  smaller  than  50  cps  to  better  define  the  spectrum  for  shallow  shots 
which  have  long  bubble  periods.  In  all  cases,  the  integration  was  carried 
out  to  the  end  of  the  positive  phase  of  the  last  bubble  pulse  observed 

on  the  tape  records.  The  number  of  bubbles  integrated  is  given  for  each 
shot  in  ^'dole  2. 


E(f)  =  — | A( f ) |  where  A(f)  is  the  amplitude  spectrum  described  in 
referen?8  (e)  and  pc  =  I.506  x  105„ 
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(3)  The  energy  in  octave  bands  from  low  frequencies  up  through  the 
8-J.6  kcs  band  was  also  computed  on  the  II#4  7090  for  both  cases— the  chock 
wave  alone  and  with  several  pulses.  The  lower  octave  band  energies  were 
included  since  it  was  of  interest  to  determine  the  rate  at  which  the 
energy  decreased  with  decreasing  frequency  for  frequencies  less  than  the 
bubble  period  frequency. 

3*3  Accuracy  of  DTMB  Computer  Data  Format  Translator.  Since  the  CTFT 
had  not  been  used  previously  for  analyses  of  explosion  pulses,  several 
checks  were  run.  In  one  instance,  the  same  record  was  digitized  twice 
on  different  days. 

■Hie  agreement  between  pairs  of  frequency  spectra  was  found  to  be 
quite  good  at  low  and  mid- frequencies,  except  that  even  slight  shifts  in 
the  baseline  affected  the  lowest  frequencies  noticeably.  At  the  highest 
frequenci’c,  differences  were  found  in  the  slopes  of  the  spectral  den~ity 
curves;  these  were  attributed  to  the  sampling  rate  which  was  not  fast 
enough  in  this  region. 

This  matter  is  discussed  in  more  detail  in  Appendix  A. 

.n  additics;  an  analysis  of  the  oscilloscope  data  for  a  few  shots 
was  made  in  order  to  obtain  an  overall  check  of  the  spectrum  results  digi¬ 
tized  by  the  CDFT.  Data  reduction  of  the  oscilloscope  records  was  done 
as  outlined  in  section  4  of  reference  (e).  The  two  methods  gave  good 
agreement . 


4.  ENERGY  SPECTRUM  OF  THE  TOTAL  PULSE 

4.1  Effect  of  Number  of  Bubbles  Integrated.  Figure  1  shows  two  spectra 
computed  from  the  same  digital  tape,  in  one  case  the  computation  was 
stopped  at  the  end  of  the  positive  phase  of  the  first  bubble  pulse,  and 
in  the  other  the  computation  was  stopped  after  the  third  bubble  pulse. 

There  appear  to  be  two  effects  attributable  to  the  difference  in 
integration  time.  First,  the  maximum  energy  shifts  to  higher  frequencies 
for  the  integration  through  the  larger  number  of  bubble  pulses.  Second, 
the  spectrum  computed  through  one  bubble  pulse  is  smoother  than  that 
through  three,  where  the  second  and  third  oscillations  are  somewhat  dis¬ 
torted.  Both  effects  are  due  to  the  fact  that  the  bubble  pulse  amplitudes 
and  periods  decrease  with  successive  pulses. 

4.2  Effect  of  Charge  Weight.  To  examine  the  effect  of  charge  weight  on 
the  spectrum,  a  comparison  is  made  in  figure  2  among  the  speclra  of  three 
TNT  charges  fired  at  nearly  4400  feet  and  weighing  1,  8,  and  57  pounds. 

The  1-  and  57-pound  shots  were  integrated  through  three  bubbles  and  the 
8-lb  shot  was  integrated  through  two  bubbles. 

1  \ 

In  Figure  2  the  energy  level  has  been  reduced  by  the  factor  (weight  ^ 
after  Weston  (reference  a),  and  the  frequency  has  been  reduced  by  multiplying 
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a 


by  (weight  5).  It  is  observed  that  these  reduction  factors  result  in 
good  agreement  among  spectra  of  varying  charge  weight  ranging  from  1  to 
57  pounds  where  the  charge  composition  and  burst  depth  are  the  same. 

^•3  Effect  of  Range-Depth.  Since  the  shots  were  recorded  near  the  sur¬ 
face,  the  effects  of  burst  depth  and  propagation  distance  upon  the  spectra 
are  difficult  to  separate  (reference  b).  The  effect  of  depth-range  upon 
the  spectrum  is  illustrated  in  Figure  3  in  which  a  plot  of  the  spectra 
of  three  1-lb  HBX-3  charges  where  the  vertical  charge-to-gage  range,  R  , 
is  equal  to  115^  feet,  U335  feet,  and  96j7  feet  respectively.  v 

In  Figure  3  the  energy  level  has  been  reduced  by  R  “,  a  spherical 
spreading  factor.  The  spectra  are  affected  by  depth,  as  expected.  The 
low  frequency  differences  are  caused  by  the  effect  of  depth  on  the  bubble 
pe.Aod;  the  maximum  energy  is  at  a  frequency  proportional  to  z  5/»f  where 
Zq  is  the  hydrostatic  depth.  The  apparent  attenuation  of  higher  frequencies 
with  range  may  or  may  not  be  real;  resolution  of  this  question  requires 
further  analysis  of  these  data. 

Effect  of  Charge  Composition.  In  Figure  U,  three  one-lb  charges  of 
different  compositions  f^red  at  a  vertical  range  of  about  1150  feet  are 
compared.  The  relatively  flat  slope  of  -U.5  dB  for  HBX-3  in  the  7-l6  kes 
band  is  probably  an  artifact  caused  by  the  sampling  rate  discussed 
previously  (Section  3«3  and  Appendix  A ).  The  slopes  of  the  other  two 
spectra  are  -9  dB  arvl  -10  dB  per  octave  for  TNT  ard  pentolite,  respectively. 
The  TOT  and  pentolite  spectra  appear  to  agree  rather  well.  In  contrast 
the  HBX-3  spectrum  has  two  differences: 

(1)  The  naxitnum  peak  and  corresponding  peaks  have  shifted  to  lower 
fre<-  iencies  because  of  the  longer  first  bubble  period. 

(2)  The  amplitude  of  oscillation  of  the  spectrum  for  HBX-3  is  not 
so  great  as  for  TOT  or  pentolite.  This  is  attributed  to  bubble  pulse 
pressures  damping  out  at  a  faster  rate  than  those  of  TOT  and  pentolite. 

The  slopes  of  the  three  spectra  are  similar  at  frequencies  below 
the  bubble  period  frequency. 

Uh 

Figure  5  shows  the  maximum  energy  spectral  density  reduced  by  W  ' 
plotted  versus  vertical  range  for  all  the  1965  data.  The  scatter  in  the 
data  Is  such  that  no  statistically  significant  conclusions  can  be  drawn. 
However,  the  HBX-3  points  would  average  perhaps  2  dB  higher  than  the  TOT 
points.  The  TOT  points  average  about  2  dB  lower  than  those  repjrted  by 
Christian  and  Blaik  (reference  b)  for  TNT  and  pentolite  combined.  Finally, 
the  NItramex  data  appear  to  be  the  lowest. 


* 

\ 


i 


5.  OCTAVE  BAND  ANALYSIS 

Using  Simpson's  Rule,  the  energy  in  ten  octave  bands  starting  with 
the  15-30  cps  band  and  including  the  8-l6  kes  band  was  computed  on  the 
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IBM  7090  computer  for  all  the  1965  shots.  This  analysis  was  performed 
for  the  positive  phase  of  the  shock  wave  (except  for  the  lowest  bands) 
and  also  for  the  pulse  including  several  bubbles,  and  is  presented  in 
Tables  *  t*rvd  4,  respectively. 


6*1  Octave  Band  Energy  for  the  Shock  Wave.  Figure  6  is  a  plot  of  octave 
band  energy  of  t-'ie  positive  phase  of  the  shock  wave  for  four  HBX-3  charges, 
weighing  from  1  to  50  pounds  and  fired  at  a  vertical  range  of  about  4200 
feet.  The  experimental  energy  level  was  redviced  by  multiplying  by  the 
weight  factor,  and  multiplied  by  the  spherical  spreading  term 

relative  to  100  yards .  fR(yds )/l00j^.  Also,  the  frequency  was  multiplied 
by  the  factor  (weight1-  ■*).  This  results  ir.  the  energy  spectrum  level  for 
a  1-lb  charge  at  100  yards  range. 


The  next  four  plots  (Figures  7-10)  present  the  shock  wave  octave 
band  energy  of  all  the  1965  data,  for  TUT,  HBX-3,  pento.Ute,  and  Nitramex, 
respectively.  The  octave  band  energies  and  the  frequency  were  reduced 
by  the  same  factors  as  in  Figure  6,  and  each  curve  was  obtained  by  drawing 
by  eye  an  average  line  through  the  data  for  each  of  the  nominal  ranges. 

In  general,  the  scatter  about  each  of  the  curves  drawn  is  about  the  same 
aj  that  shown  in  Figure  6.  However,  for  certain  depths*  the  relatively 
few  data  points  and  the  scatter  of  the  d«tn  in  the  8-l6  kes  band  resulted 
In  eurv?-..  which  are  not  reliable  ir.  the  high  frequency  region. 


At  the  very  low  frequency  end  of  the  spectrum,  the  energy  level 
approaches  the  shock  wave  impulse.  In  this  region,  the  energy  decreases 
with  increasing  depth  and  this  is  characteristic  of  all  data.  This  is 
expected  since  the  shock  wave  impulse  decreases  faster  than  range"1 
(reference  f)  for  all  compositions.  The  spread  in  the  reduced  spoctra 
decreases  with  increasing  frequency  until  a  minimum  is  obtained  at  about 
2000  cps  x  lbs1'’.  As  the  frequency  increases  further,  the  spread  in 
energies  increases  again;  the  energy  in  general  decreases  with  increasing 
depth  as  occurred  at  the  low  frequencies. 

The  change  in  the  shock  wave  spectra  at  the  low  frequency  end  as 
the  depth  is  varied  is  a  consequence  of  the  change  in  the  waveform  with 
depth.  As  the  depth  increases,  the  shock  wave  duration  decreases,  resulting 
in  a  shorter  pulse  vdiich  will  produce  less  energy  in  the  lower  frequencies. 

The  HBX-3  spectra  show  1  to  3  dB  higher  energy  than  TJfT  at  the  lowest 
frequencies  for  all  ranges.  Since  the  spectrum  approaches  the  impulse 
at  zero  frequency,  this  higher  energy  level  for  HBX-3  is  caused  by  higher 
shock  wave  impulse  (reierence  f).  In  general,  tne  TNT  and  HBX-3  data  arc 
in  good  agreement  at  all  but  the  lowest  frequencies.  Pentolite  is  1  dB 
to  3  dB  higher  than  TlfT  for  all  frequencies  at  the  three  ranges  for  which 
data  exist.  Nitramex  is  2  dB  or  3  dB  lower  than  TNT  at  most  frequencies 
for  corresponding  depths.  Nitramex  has  a  slope  of  -6  dB  per  octave  in 
the  high  frequencies;  this  is  a  relatively  high  slope. 


*  TNT  data  at  nominal  depths  of  1200  ft  and  4500  ft;  HBX-3  data  at 
nominal  depths  of  1200,  2000,  4^00,  and  10,000  ft;  Nitramex  data  at 
the  nominal  depth  cf  3000  ft. 
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There  are  two  slight  anomalies  in  the  HBX-3  data:  the  2000  ft  data 
are  consistently  1  or  2  dB  higher  than  the  1200  ft  data,  and  the  energies 
for  the  one  3000  ft  depth  shot  fired  are  higher  than  the  1200  ft  depth 
energies.  Again  these  comparisons  are  based  on  statistically  inadequate 
data. 

5*2  Octave  Band  Energy  for  the  Total  Pulse.  The  next  four  plots  (Figures 
11- 14)  show  the  octave  band  energy  for  the  total  pulse  (shock  wave  plus 
one  or  more  bubble  pulses  as  indicated  in  Table  2)  far  the  four  comnosi- 
tions.  Again  la*  energy  level  is  reduced  by  veightV3  and  multiplied  by 
the  spherical  spreading  factor;  the  frequency  is  reduced  by  welght^/3  as 
in  Figures  7~10.  The  octave  band  energy  density  points  were  connected 
by  straight  lines  and  the  maximum  value  is  an  extrapolation  of  these 
curves  at  the  bubble  period  frequency. 

It  is  observed  in  the  four  sets  of  data  that  there  is  no  systematic 
change  in  the  slope  of  the  line  for  frequencies  less  them  the  bubble 
period  frequency.  An  accurate  determination  of  the  slopes  for  these 
frequencies  could  not  be  made  since  in  general  there  were  only  two  octave 
bands  below  the  bubble  period  frequency.  Furthermore ,  the  difficulty 
in  determining  the  value  of  the  baseline  (discussed  in  Section  3»l)  may 
have  resulted  in  large  variations  of  energy  for  the  lowest  bands.  However, 
the  curve  in  this  region  usually  has  a  steeper  slope  (8  dB  or  9  dB  per 
octave)  than  the  slope  of  3  dB  or  4  dB  per  octave  reported  by  Christian 
and  Bluik  (reference  b),  and  is  believed  to  be  more  realistic,  since  it 
is  based  on  considerably  more  data  and  a  more  closely  spaced  frequency 
analysis. 


In  comparing  the  data  for  different  charge  compositions  at  the  same 
ranges,  it  was  observed  that  TNT  and  HBX-.3  have  comparable  energy  levels 
for  reduced  frequencies  greater  than  1  kc'W1^.  However  in  the  low  fre¬ 
quency  region,  the  maximum  energy  for  KBX-1  is  greater  than  that  for  TNT 
and  occurs  at  lower  frequencies.  Pentolite  has  a  2  dB  to  4  dB  higher 
energy  level  than  TlfT  for  most  frequencies.  It  must  be  noted  that  there 
is  only  one  pentolite  shot  and  only  one  weight— 1  pound — at  each  of  the 
three  ranges.  Nitramex  was  observed  to  have  a  2  dB  to  4  dB  lower  energy 
level  than  TffT  for  all  frequencies;  this,  too  in  based  on  scanty  data. 


6.  CONCLUSIONS 

(1)  The  good  agreement  previously  found  between  the  reduced  total 
prise  energy  _pectra  of  1-  and  10-pound  charges  has  been  shown  to  hold 
up  to  57-p«  id  charges. 

(2)  The  total  pulse  energy  spectra  cf  charges  fired  at  various 
vertical  renges  between  500  and  14,000  feet  vary  with  depth  as  expected 
from  previous  work. 

(3)  The  total  pulse  energy  spectra  of  TNT  and  pentolite  are  almost 
identical;  the  HBX-3  spectrum  show*:  a  maximum  (and  succeeding  peaks)  at  a 
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lower  requenry,  less  oscillation,  and  smaller  attenuation  at  the  highest 
frequencies.  The  energy  of  the  maximum  was  not  statistically  different 
for  the  four  compositions  fired;  however,  HBX-3  gave  the  highest  values 
ana  Ifitramex  the  lowest. 

(M  The  octave-band  energy  for  the  shock  wave  alone  and  for  the 
total  pulse  show  the  effect  of  depth  in  sharpening  the  shock  waveform  and 
decreasing  the  bubble  period.  The  low  frequency  drop  off  of  8-9  dB  per 
octave  for  the  total  pulse  energy  is  believed  to  be  a  better  value  than 
the  previous  one. 
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TABLE  1 

COMPARISON  OP  PREVIOUS  AND  CURRENT  EXPERIMENTAL  CONDITIONS 


Christian  and  Blaik 
(reference  b) 

February  19^5 

Sea  Trials 

Charge  Weight 

0.4  to  10  lbs 

1  to  88  lbs 

Charge  Depth 

7000-22,000  ft 

500-14,000  ft 

Charge  Composition 

TNT,  Pentolite 

TNT,  Pentolite,  HBX-3, 
Nitraaex 

Position  of  Recording 
Gages 

Vertically  above 
charge  near  surface 

Vertically  above 
charge  near  surface 

Pulse  Sampling  Method 

Pressures  sampled  at 
discrete  time  intervals 
from  paper  and  film 
records 

Pressures  from  analog 
tape  data  electronically 
sampled  and  digitized 
using  Computer  Data 
Format  Translator 

Length  of  Pulse 

Analyzed 

To  minitnjtm  pressure 
or  end  of  negative 
phase  after  second 
bubble  pulse 

To  end  of  positive 
phase  of  last 
observable  bunoie 
pulse 

Lowest  Octave  Band 

250-500  cps 

15-30  cps 
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TABLE  2 


SHOT  STATISTICS 


Charge 

Weight 

(lb) 

Charge 

Composition 

Nominal 

Depth 

(ft) 

Burst 

Depth 

(ft) 

Shot 

Number 

Number  of  Bubble 
Pulses  Included  In 
Spectrum  Computation 

1 

TNT 

500 

517 

22 

1 

500 

575 

48 

2 

500 

607 

33 

1 

800 

873 

49 

2 

800 

893 

67 

2 

1,200 

1,158 

42 

3 

1,200 

1,241 

70 

2 

2,000 

1,983 

27 

3 

3,000 

2,963 

45 

3 

3,000 

3,032* 

10 

3 

4,500 

4,411 

24 

3 

4,500 

4,502* 

52 

3 

7,000 

7,402* 

77 

3 

1 

PENTOLITE 

500 

552 

64 

3 

1,200 

1,377 

68 

2 

7,000 

^,621* 

5 

3 

1 

HBX-3 

1,200 

1,331 

46 

3 

2,000 

1,891 

65 

4 

4,500 

4,228* 

51 

3 

4,500 

4,517 

56 

3 

10,000 

9,824 

73 

4 

8 

TNT 

500 

430 

26 

1 

500 

571 

31 

1 

800 

892 

12 

1 

1,200 

1,139 

72 

3 

1,200 

1,237 

43 

2 

1,1300 

1,250 

69 

2 

2,000 

1,871 

21 

3 

2,000 

2,031 

35 

3 

2,000 

2,051* 

4 

2 

3,000 

2,924 

47 

3 

4,500 

4,305 

32 

3 

4,500 

4,372 

25 

2 

4,500 

4,542 

55 

3 

14,000 

13,300* 

18 

3 

14,000 

13,400 

36 

3 

8 

HBX-3 

1,200 

1,182 

44 

3 

1,200 

1>  358 

11 

2 

4,500 

4,202 

50 

2 

10,000 

9,849 

71 

3 

14,000 

13,540* 

75 

3 

*  Depth  determined  from  bubble  period. 
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TABLE  2 


Charge 

Weight 

Charge 

Composition 

SHOT  STATISTICS  (cont’d) 

Nominal  Burst  Shot 

Depth  Depth  Number 

_ £*i _ m _ 

Number  of  Bubble 
Pulses  Included  In 
Spcctrutt  Conroutation 

57 

TNT 

2,000 

1,844 

8 

2 

2,000 

1,923 

29 

2 

4,500 

4,282 

37 

2 

4,500 

4,413 

20 

3 

7,000 

6,608* 

3 

3 

7,000 

7,147* 

79 

3 

14,000 

14,340 

62 

3 

50 

HBX-3 

2,000 

1,313 

63 

2 

3,000 

2,891 

78 

3 

4,500 

4,443 

53 

2 

7,000 

7,498 

57 

3 

10,000 

9,547 

59 

2 

88 

NITRAMEX 

3,000 

2,793 

38 

2 

3,000 

2,900 

19 

1 

7,000 

7,051 

39 

3 

10,000 

10,308 

40 

3 
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SPECTRA  Ox’  RJLSE  INCLUDING  SEVERAL  BUBBLES  (continued) 
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SHOT  5  1-LB  PENTOUTE  BURST  DEPTH  =6620  FT. 
(DATA  OM  SAME  DIGITAL  TAPE  FILE) 
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FiG.  1  EFFECT  OF  NUMBER  OF  BUBBLES  INTEGRATED  UPON 
ENERGY  SPECTRUM 


TNT  CHARGES  FIRED  AT  NEARLY  4400  FT. 
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FIG.  2  EFFECT  OF  CHARGE  WEIGHT  ON  ENERGY  SPECTRUM 


-LB  HBX-3  CHARGES 
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FIG.  3  EFFECT  OF  RANGE -DEPTH  UPON  ENERGY  SPECTRUM 
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FIG.  4  EFFECT  OF  CHARGE  COMPOSITION  UPON  ENERGY  SPECTRUM 
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FIG.  5  TOTAL  PULSE  MAXIMUM  SPECTRAL  ENERGY  DENSITY 
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FIG.  6  SHOCK  WAVE  OCTAVE  BAND  ENERGY  FOR  HBX-3  CHAR. 


3MINAL  BURST  DEPTH 


NOLTR  66-  128 


o' 


n 

\ 


i 

.  IS) 
>-  Q- 

U  U 
Z 

UJ  to 


0  - 


cs. 

u. 


O  7T. 


H 

0 


u 

D 

Q 

LU 

n C 


£ 

X 

u 

z 

UJ 

0 

UJ 

O' 

u. 


L 

o 


- 1 - 1 - ! _ L__ 

°  ©  o  O 

CM  CO 

1  I  I 

(SGA  001  iV  308VHD  01- t  803 

(‘j  <0/3/0/ 083  l  38  8P  )  '  (  -22L\  C/frM 

2  \(13 )  8/  X  3 


J 


§ 


o 


28 


FIG.  7  SHOCK  WAVE  OCTAVE  BAND  ENFRGY  FOR  TNT  CHARGES 


+  10  r  NOMINAL  BURST  DEPTH 


NCLTR  66-128 


29 


FIG.  8  SHOCK  WAVE  OCTAVE  BAND  ENERGY  FOR  HBX-3  CHARGES 


NOMINAL  BURST  DEPTH 
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FIG.  9  SHOCK  WAVE  OCTAVE  BAND  ENERGY  FOR  PENTOLITE  CHARGES 


NOMINAL  BURST  DEPTH 
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FIG.  10  SHOCK  WAVE  OCTAVE  BAND  ENERGY  FOR  NIT.1AMEX  CH,^GES 


NOMINAL  BUST  DEPTH 
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FIG.  11  TOTAL  PULSE  OCTAVE  BAND  ENERGY  FOR  TN"  CHARGES 
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FIG.  12  TOTAL  PULSE  OCTAVE  BAND  ENERGY  FOR  HBX-3  CHARGES 


NOMINAL  BURST  DEPTH 


NOLTR  66-128 


( SQA  001  iV  39HVHD  81-1  SOd 


SdD/zWD/9J)3l  3d  8P )  x 


34 


FIG.  13  TOTAL  PULSE  OCTAVE  BAND  ENERGY  FOR  PENTOLITE  CHARGES 
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APPFNDIX  A 

REPRODUCIBILITY  OF  SPECTRA.  FROM  DIGITIZED  DATA 


In  order  to  test  the  ETMB  Computer  Data  Format  Translator  with  our 
analog  tapes,  repeat  runs  of  digitization  were  made  and  the  results  com¬ 
pared.  This  was  done  for  a  record  which  was  analyzed  only  for  the  positive 
shock  wave  pulse,  and  for  another  record  analyzed  through  the  first  bubble 
pulse. 


Figure  A-l  illustrates  typical  energy  spectra  of  the  positive  phase 
of  the  shock  wave.  The  two  plots  shown  are  of  the  same  shot,  No.  US,  a 
1-lb  lfBX—  ^  charge  fired  at  1330  feet,  the  difference  being  that  the  analog 
data  was  digitized  on  different  days.  Tliis  was  done  to  check  the  repro¬ 
ducibility  of  the  iTTMB  equipment.  The  spectrum  becomes  level  at  the  lower 
frequencies  and  approaches  the  impulse*  at  zero  frequency.  The  reproduci¬ 
bility  of  the  shock  wave  spectra  at  low  and  raid  frequencies  is  observed 
to  be  quite  good. 

The  spectrum  of  an  exponential  pulse  should  ideally  have  a  -6  dB 
per  octave  slope  in  the  high  frequencies  (reference  g).  However,  in 
Figure  A-l  the  slopes  of  the  spectra  are  -3  dB  per  octave  and  -11  dB  per 
octave,  respectively.  Although  the  sampling;  was  at  33  usee  intervals  in 
both  cases,  the  sampling  was  probably  initiated  at  different  times, 
resulting  in  digitizing  different  portions  of  the  pressure-time  data. 

In  this  particular  case,  where  th'  ...mes  are  the  order  of  50  usee, 

P**  about  150  usee,  and  the  sampling  interval  33*3  usee,  reading  differ¬ 
ent  pressures  on  the  rise  and  near  the  peak  could  result  in  different 
slopes  in  the  spectrum  at  the  highest  frequencies.  This  was  verified 
by  reading  the  same  pressure-time  Visicorder  record  on  the  Telereadex  in 
such'  a  way  that  the  sampling  was  controlled  by  starting  the  sampling  at 
different  times  relative  to  the  discontinuous  rise  of  the  shock  wave.  It 
becomes  apparent  that  the  sampling  should  have  been  at  closer  intervals; 
however,  the  available  lb. 6  usee  sampling  rate  was  not  used  because  it 
introduced  too  much  loise,  as  mentioned  before. 

Figure  A -2  shows  e.-amples  of  energy  spectra  of  a  pulse  integrated 
to  the  end  of  the  positive  phase  of  the  first  bubble  pulse.  Here  again 
the  same  analog  data  was  digitized  twice.  Ar  expected  for  oscillating 
functions,  the  maximum  energy  occurs  at  the  oscillating  frequency;  i.e., 
at  the  bubble  period  frequency  ( 35  cps )  and  then  continues  to  oscillate 
at  integral  multiples  of  this  frequency.  At  500  cps,  the  interval  in 
which  the  spectrum  was  computed  was  changed  from  5  to  50  cps;  therefore, 

*  Impulse  is  defined  as  JQ  p(t)dt,  where  T  is  the  duration  of  the 
pressure  pulse,  p(t),  being  integrated. 

**  q,  the  time  constant  of  the  pressure  pulse,  is  defined  as  the  time 
where  the  pressure  falls  to  l/e  of  its  maximum  value,  p  . 
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the  35  cps  oscillations  are  not  veil  defined  for  frequencies  greater  than 
300  cps.  It  is  apparent  that  the  spectra  are  in  gocd  agreement  from  the 
bubble  fundamental  frequency  to  about  6-8  kcs.  In  the  high  frequencies, 
the  slight  differences  in  slope  are  again  attributed  to  the  relatively 
coarse  sampling  rate  and  the  uncontrolled  starting  point  of  the  sampling. 

Although  the  characteristics  of  the  tvo  curves  fron  about  20  cycles 
to  about  5  cycles  are  similar,  the  discrepancies  are  3  to  5  dB-  This 
large  discrepancy  is  probably  due  to  the  different  baselines  calculated, 
since  a  slight  shift  in  the  baseline  could  produce  relatively  large 
differences  in  the  Impulse. 
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FIG.  A-l  REPRODUCIBILITY  OF  SHOCK  WAVE  SPECTRUM 


SHOT  12  -  8- LBS  TNT  BURST  DEPTH  =  890  FT 


NOLTR  66-128 


I _ l _ i _ i _ i 

o  o  o  o  o 

CM  —  —CM 

+  +  II 

(S<D/2WD/983  l  38  SP) 
AKSN3Q  1V81D3dS  V38V  ilNn/A983N3 


A-4 


Unclassified 


DOCUMENT  CONTROL  DATA  •  R  &  D 

.'>«•«  unis  r/ds  (n  ation  t >/  title,  bods  vt  nb^tntet  and  indexing  annotation  must  be  entered  when  the  overall  report  is  t  lassihed) 

'  originating  activity  (Corpora te  a  thor) 

2a.  report  security  Classify  ation 

U.  S.  Naval  Ordnance  Laboratory 

Unclassified 

White  Oak,  Silver  Spring,  Maryland 

3  REPORT  TITLE 

Spectrum  and  Octave  Band  Analysis  of  Pressure  Pulses  From  Deep  Underwater 

Explosions 

a  descriptive  noi  £S  fTVpe  0/  repoK  And  inr/usive  dslesi 

S  aij  ThoRiS)  (First  name,  middle  initial,  last  name) 

Mary  Alice  Genau 

6  REPOR  T  O  *  T  E 


Bm  CONTRACT  OR  GRANT  NO 

b  DflOj£C  T  no  rpask  NOL-785/ARPA 

c. 

d. 


5  October  19 66 


10  DISTRIBUTION  STATEMENT 


7 a.  TOTAL  NO  0»  PAGES 

50 


90  ORIGINATOR’S  REPORT  NUMBERlS) 


NOI/TR  66-128 


9fc.  OTHER  REPORT  NO  I  Si  (Any  other  number*  (hot  may  bo  assigned 
thi •  report ) 


u. supplementary  notes 


13  ABSTRACT 


Distribution  of  this  document  is  unlimited. 


12  SPONSOHINC  MILITARY  ACTIVITY 


Advanced  Research  Project  Agency 


Fourier  spectra  were  computed  on  the  IHM  7090  for  analog  tape 
recorded  pressure  pulses  of  underwater  explosions  fired  at  sea  in 
February  1965.  Depths  ranged  from  500  to  14,000  feet;  charges  weighed 
1  to  88  pounds;  the  compositions  fired  were  TOT,  pentolite,  HBX-3,  and 
Nitramex.  Reduced  spectra  of  charges  weighing  up  to  57  pounds  agreed 
with  previous  results  from  1  and  10  pound  charges  at  the  same  depths 
Only  slight  differences  due  to  composition  were  found. 


DD  ,Fr\,1473 


(PAGE  1 


S/N  0101-807. sect 


Securtv  Classification 


Urr’  -water  ’"xplosion 

She  Waves 

Frequency  Sp<  -tram 

’Underwater  Sound 

Bubble  Pulses 

TUT 

H3X-3 

Pentollte 

Nltranex 


DD  ."“.1473 

(PAGE  2) 


( BACK  ) 


